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Abstract 

The chiral helimagnet Cri/3NbS2 has been investigated by magnetic, transport and thermal 
properties measurements on single crystals and by first principles electronic structure calculations. 
From the measured field and temperature dependence of the magnetization for fields applied per- 
pendicular to the c axis, the magnetic phase diagram has been constructed in the vicinity of the 
phase transitions. A transition from a paramagnetic to a magnetically ordered phase occurs near 
120 K. With increasing magnetic field and at temperatures below 120 K, this material undergoes 
transitions from a helimagnetic to a soliton-lattice phase near 900 Oe, and then to a ferromagnetic 
phase near 1300 Oe. The transitions are found to strongly affect the electrical transport. The resis- 
tivity decreases sharply upon cooling near 120 K, and the spin reorientation from the helimagnetic 
ground state to the commensurate ferromagnetic state is evident in the magnetoresistance. At 
high fields a large magnetoresistance (55 % at 140 kOe) is observed near the magnetic transition 
temperature. Heat capacity and electronic structure calculations show the density of states at the 
Fermi level is low in the magnetically ordered state. Effects of spin fluctuations are likely important 
in understanding the behavior of Cr 1 / 3 NbS2 near and above the magnetic ordering transitions. 



PACS numbers: 75.30.Kz, 75.10.-b 



I. INTRODUCTION 



Complex magnetic textures have recently been observed in some magnetic materials crys- 
tallizing in the cubic space group P2 1 3P s MnSL Feo.5Coo.5Si, FeGe, and Cu 2 OSe0 3 are those 
materials where the ground state magnetic structure is found to have a modulated helimag- 
netic arrangement of spins incommensurate with the underlying crystal lattice. In addition 
to this, another remarkable similarity in these materials is that in a certain range of magnetic 
fields and at temperatures just below the transition temperature, a distinct skyrmion lattice 
phase exists. A skyrmion is a magnetic vortex, a structure reminiscent of the Abrikosov vor- 
tices in type II superconductors.™^ It is to be noted that these similarities have been observed 
irrespective of the nature of the electrical conductivity of the materials. MnSi and FeGe are 
metals, Feo.5Coo.5Si is a small band gap semiconductor and Cu20Se03 is an insulator. The 
first three have ferromagnetic ordering and the latter one has ferrimagnetic ordering in a 
certain plane that modulates along a particular crystallographic direction to form a helical 
ground state. The helix eventually gets destabilized to long ranged ferromagnetic and ferri- 
magnetic ordering respectively with the application of magnetic field. A common property 
all the materials share is the symmetry of the crystal structure. The remarkable feature of 
space group P2i3 is that it lacks an inversion center. The ground state chiral helimagnetic 
ordering in the non-centrosymmetric space group has been well understood to be due to the 
competition between the symmetric exchange interaction that favors a parallel arrangement 
of spins and the antisymmetric Dzyaloshinsky-Moriya(DM) interaction that favors a perpen- 
dicular arrangement of spins, which results in a modulated helimagnetic spin structure.^^ 
The origin of the skyrmion phase, however, is still unclear on a general level. Theoretically, it 
has been predicted to occur in a wide range of non-centrosymmetric mag 

nets EHll5] However, 

experimentally they have been observed in materials belonging to only one space group. At 
the same time, recent studies have shown that the non-centrosymmetric magnets provide a 
platform for the observation of a magnetic blue phase, a novel magnetic state, formed by 
the twisting of the chiral helices, reminiscent of the blue phase observed in nematic liquid 
crystals.™ 3 

Cr!/ 3 NbS2 is a material that has all the common properties of the materials mentioned 
above except that it bears a different crystallographic symmetry. It crystallizes in a non- 
centrosymmetric hexagonal space group P6322 belonging to the point group Dgp3 one of 
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the symmetry groups in which ferromagnets have been theoretically predicted to have a 
skyrmion lattice™ It orders magnetically into a helimagnetic ground state with a period 
of about 480 A below about 120 - 127 K and is known to have metallic behavior below 
room temperature.^^ The spins are arranged ferromagnetically in the afr-plane and the 
helix is along the c axis. Effect of magnetic field in the a6-plane has been found to be 
dramatic with a metamagnetic transition observed at an applied field near 1200-1500 Oe.^ 
A recent study conducted on a thin crystal with small angle electron diffraction and Lorentz 
force microscopy has shown that a magnetic field applied perpendicular to the direction of 
the helix destabilizes the helical structure gradually into a soliton lattice, a nonlinear pe- 
riodic magnetic state, with an eventual incommensurate-to-commensurate transition into a 
ferromagnetic state at the critical field of 2300 Oe.^ Manipulation of the spin spiral with 
magnetic field has generated interest in this material for spintronics applications.^^ 

The first helimagnetic ground state in this material was experimentally observed by 
Miyadai et alW^ by means of small angle neutron scattering conducted on a powder sample, 
which was in accordance with the prediction of the helimagnetic structure by Moriya and 
MiyadaP^ based on the magnetic measurements and the subsequent theoretical interpreta- 
tion in terms of the DM interaction in a non-centrosymmetric magnet. In earlier studies 
Parkin et 6tZJ 22 * 23 * found it to order ferromagnetically below 120 K with the basal plane being 
the easy axis. 

In this paper we present experimental results obtained on high quality single crystals of 
Cr!/ 3 NbS2 together with results from electronic structure calculations. Temperature depen- 
dence of DC susceptibility and the magnetization curves are measured both parallel and 
perpendicular to the c axis. Likewise, the temperature dependence of electrical resistivity, 
specific heat capacity, thermal conductivity and Seebeck coefficient in zero magnetic field 
and their behavior in magnetic fields applied perpendicular to the c axis are investigated. 
Results of band structure and density of states (DOS) calculations are presented both in the 
non-magnetic and magnetic states. 

II. CRYSTAL CHEMISTRY 

Cri/ 3 NbS2 crystallizes in the NbsCoSg, hp20 structure type™It forms a hexagonal layered 
structure as depicted in Figure pTa) containing of 20 atoms per unit cell. Twelve sulfur atoms 
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FIG. 1: (Color online) a) Hexagonal crystal structure of Cr 1/ / 3 NbS2- The intercalated chromium 
atoms occupy the octahedral interstitial holes between the trigonal prismatic layers of 2H-NbS2. 
b) Structure of Cr 1 / 3 NbS2 emphasizing the layers of Cr atoms in the afr-plane. The Cr-Cr distance 
is shortest in the plane. The dotted lines show the unit cell. The biggest (orange) balls represent 
Cr atoms, the medium (green) balls represent Nb atoms and the smallest (black) balls represent S 
atoms. Nb atoms are in two inequivalent sites labeled as Nbl and Nb2. 

occupy the general site. There are six niobium atoms in two inequivalent positions - 4f and 
2a. The Cr atoms are intercalated in the octahedral holes (2c sites) between the trigonal 
prismatic layers of 2H-NbS2- The unit cell parameters obtained from the x-ray data collected 
from the powdered crystals are a = 5.741 A and c = 12.101 A which are in good agreement 
with the reported values. EED The Cr-Cr distance is closest within the layer in the afr-plane 
(5.741 A). Along the c axis, the nearest Cr-Cr distance is 6.847 A. Cr^Nb^ belongs to a 
large family of layered materials that are formed by intercalation of transition metal elements 
within the layers of the first row transition metal dichalcogenides. The structure depends 
on the amount of intercalation in the layer. In the absence of intercalated atoms the host 
transition metal dichalcogenide layers are coupled through weak van der Waals bond. The 
intercalation of a transition metal strengthens the bonding with the possibility of charge 
transfer from the intercalated atoms to the transition metal atoms in the layer which brings 
about strong changes in the electronic structure and corresponding changes in the electrical 
transport and magnetic properties.®' 
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III. EXPERIMENTAL DETAILS 



Crystals of Cr 1 / 3 NbS2 were grown by chemical vapor transport using iodine as the trans- 
port agentP^. Single crystals of various sizes with dimensions as large as 8 mm x 7 mm x 
1 mm were obtained. X-ray diffraction from powdered crystals confirmed the P6322 struc- 
ture and revealed no impurity phases. Energy dispersive x-ray spectrometer (EDS) results 
along with the magnetic transition temperature clearly distinguish these crystals as having 
1/3 intercalation of Cr atoms.^^ No iodine was detected in the crystals by the EDS mea- 
surements. Magnetic properties were measured using a Quantum Design magnetic property 
measurement system (MPMS). Resistivity, magnetoresistance, thermal conductivity and 
Seebeck coefficient were all measured in a Quantum Design physical property measurement 
system (PPMS). A four wire configuration with platinum wires and Epotek H20E silver 
epoxy were used for the the resistivity and magnetoresistance measurements. 



IV. RESULTS AND DISCUSSION 



A. Magnetic properties 



Figure [2] shows M/H as a function of temperature with the applied magnetic field per- 
pendicular [Fig. [2^a)] and parallel [Fig. [2^b)] to the c axis. With the magnetic field applied 
perpendicular to the c axis, at the lower field (1 kOe), a small kink is observed at about 117 
K. This kink is not observed at 10 kOe. No such kink is observed in M/H vs. T when the 
magnetic field is applied parallel to the c axis. These results are consistent with previous 
measurements by Miyadai et al. 21 As the material is known to have helimagnetic ordering 
in lower fields with the helix directed along the c axis,^ the kink may represent the onset 
of the helical state. Such a kink has also been observed in other known helimagnets MnSp^ 
FeGeP and Fei- X Co^SP^. Cri/ 3 NbS2 follows Curie- Weiss behavior at higher temperatures. 
The inset in Fig. 2J^a) shows the Curie- Weiss fit oi x 1 — qSe — ^° ^ ne high temperature part 
of data measured with the applied field of 1 kOe. The parameters obtained are the Curie 
constant C = 2.4 K cm 3 mol _1 F.U. -1 and the Curie- Weiss temperature 9cw — 127 K. The 
effective moment per mole of chromium atoms, p e ff, calculated from the Curie constant is 
4.4 fiB, which is consistent with the values reported by Parkin et alW^ and is close to the 
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FIG. 2: (Color online) M/H as a function of temperature at the fields indicated with the magnetic 
field applied a) perpendicular to the c axis, and b) parallel to the c axis. Inset in (a) shows the fit 
to the Curie- Weiss law for the data taken at 1 kOe above 160 K. 

spin only value of 3.87 hb for Cr 3+ . 

The anisotropic nature of the low temperature magnetism in Cri/3NbS2 is evident in Fig. 
[3] Applying a field along the c axis [Fig. |3^b)] simply rotates the ordered moments out 
of the a6-plane. This results in a nearly linear increase in M with H up to H = 20 kOe. 
When the field is applied in the afr-plane [Fig. |3](a)], more complex behavior is observed, 
which will be addressed further in the following discussion. The moment saturates to the 
same value of about 3.2 \i% /Cr in both orientations, but this occurs at a much lower field 
when H is applied in the afr-plane, indicating a strong preference for the moments to remain 
perpendicular to the c axis. 

Figure [4] shows various magnetic properties near the magnetic ordering temperature with 
the magnetic field applied perpendicular to the c axis i.e. perpendicular to the helical axis. 
M/H as a function of temperature at the indicated magnetic fields are depicted in Fig. [3|a). 
It is observed that the kink due to the onset of helimagnetic ordering vanishes above 1100 
Oe. At 119 K, M vs. H increases linearly at low field. At around 500 Oe, it shows a change 
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FIG. 3: M vs. H measured at a temperature of 2 K with the magnetic field applied a) perpen- 
dicular, and b) parallel to the c axis. Note the different scales on the x-axis. 

in slope. Another slope change is observed around 1000 Oe. The region between 500 - 1000 
Oe shows a steeper slope. With the decrease in temperature, the lower end of the steeper 
slope moves towards higher field and a sharp change is observed around 1000 Oe. The mo- 
ments saturate above about 1300 Oe. These observations in the magnetic measurements are 
consistent with the helimagnetic ordering and soliton lattice phase formation in the material 
as observed by Togawa et a/P In zero applied field, below the transition temperature the 
moments in a plane perpendicular to c axis align parallel to each other while the moments 
in the plane make small angle to the parallel moments in the adjacent planes. The kink thus 
represents the effect of the small antisymmetric coupling between the planes giving rise to 
a long period helimagnetic ordering. As an external magnetic field is applied perpendicular 
to the c axis, at smaller fields (up to 500 Oe) the field is not strong enough to overcome 
the antisymmetric coupling between the planes. It can, however, turn the moments in the 
planes gradually in the direction of the field. Thus, M increases slowly with H at low fields. 
As the field gets stronger, it starts breaking the antisymmetric coupling between the planes 
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FIG. 4: (Color online) Magnetic properties of Cr^Nb^ near the transition temperature in 
low magnetic field applied perpendicular to the c axis, a) M/H as a function of temperature at 
indicated fields, b) M vs. H at selected temperatures, and c) Arrott plots. The solid lines are the 
extrapolation of the M 2 vs. M/H data at higher fields. 

and the moments in a certain plane start to rotate in the direction of the applied field. This 
results in a rapid increase of M. This region of the phase diagram is where the period of 
the helix is found to increase with the applied magnetic field forming the soliton phased A 
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large enough field (1300 Oe) aligns all the moments parallel to the applied field thus causing 
the transition to the ferromagnetic state and saturation of the moment. 

Figure [5] shows a contour plot of M vs. H and T near the phase transition region. This 
provides an useful visualization of the magnetic phase diagram for Cri/3NbS2- The white 
lines are meant to guide the eye, and follow the temperature and field dependence of the 
transition temperatures (white circles) determined from examination of individual M(H) 
curves at fixed T [see, for example, the 110 K curve in Fig. [41(b)] . At high temperature, 
Cr!/ 3 NbS2 is paramagnetic (PM). The lower set of points separates the helimagnetic (HM) 
phase from the soliton lattice (SL) phase, and is determined by the onset of the rapid increase 
in M with H. The upper set of points separates the SL phase from the ferromagnetic (FM) 
phase, and is determined by the saturation field at a particular temperature. The details of 
how the complex HM and SL phases evolve upon cooling from the PM phase have not been 
determined. Here we observed the transition from the soliton phase to the ferromagnetic 
state at field of about 1300 Oe. In the thin sample measured by Lorentz microscopy, Togawa 
et aP^ found it to occur at 2300 Oe. In our crystals, the manifestation of the metamagnetic 
transition is apparent even in the thermal and transport properties discussed below. 




T(K) 



FIG. 5: (Color online) Contours of M (H,T) determined from M vs. T measurements at fields 
from 100 Oe to 1500 Oe applied in the a6-plane. The white data points were determined from 
M vs. H measurements at temperatures from 110 to 120 K. The question mark "?" represents 
the region of transition form paramagnetic (PM) to helimagnetic (HM) phase. The details of this 
phase evolution have not yet been determined. 
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Figure |4][c) shows variation of M 2 (H, T) as a function of H/M(H } T) at different temper- 
atures in the region of the magnetic transition. These Arrott plots in the ferromagnetic state 
(at higher fields) are linear, indicating itinerant ferromagnetismP*^ These data are consis- 
tent with the result obtained from first principles calculations presented below, where the 
Stoner criterion has been found to be satisfied. The isothermal line in the Arrott plot pass- 
ing through the origin represents the transition temperature.^ Thus the Curie temperature 
(T c ) for Cr 1/3 NbS 2 is taken to be 120 K. 

B. Transport properties 

2.25 -r-5 1 




T(K) 



FIG. 6: (Color online) a) Electrical resistivity of Cr 1 / 3 NbS2 as a function of temperature mea- 
sured in the a6-plane. Inset shows the temperature derivative of the resistivity, b) Temperature 
dependence of electrical resistivity of Cr 1/ / 3 NbS2 measured in the a6-plane at indicated magnetic 
fields applied parallel to the plane. Inset shows the magnetoresistance where Ap a i, = PiiT-po- 

Electrical resistivity as a function of temperature is depicted in Figure [6|a). Interestingly, 
there is an abrupt and large change of the resistivity in the vicinity of the magnetic tran- 
sition temperature. Above Tc, up to the measured temperature of 300 K, there is a slight 
decrease in the resistivity with temperature. It contrasts with the behavior of conventional 
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metallic ferromagnets where there is a decrease in the resistivity on cooling over the entire 
temperature range and a change in slope or a kink near the magnetic transition temperature 
due to reduction in the spin disorder scattering. Hall effect measurements at 200 K with 
the current in the afr-plane and the magnetic field along the c axis gave a Hall coefficient of 
2xl0~ 2 cm 3 C -1 at 200 K (for a single band, this would correspond to 10 20 holes cm -3 ). 
This along with the measured resistivity values near 2 mQ cm above the transition, sug- 
gests the material may be described as low carrier concentration metal or heavily doped 
semiconductor. The decrease in p with increasing temperature at the highest temperatures 
investigated indicate that this is not a simple metal in the paramagnetic state. Comparison 
of the transport properties with those reported in Ref. [221 shows that the crystals studied 
here have higher resistivity and higher Hall coefficients than those previously studied. 

Figure ^b) shows the temperature dependence of the resistivity at various magnetic fields 
up to 140 kOe. As expected the effect of the magnetic field is strongest near the magnetic 
transition temperature. The sharp change observed at lower fields near Tq are observed 
to skew to higher temperatures by magnetic field, as may be expected for a ferromagnet. 
However, a strong suppression of the resistivity upon cooling is still seen at 140 kOe. Inset 
of figure |6^b) shows the magnetoresistance obtained by subtracting the resistivity measured 
at zero applied magnetic field from the resistivity measured at the magnetic field of 140 kOe. 
A very large magnetoresistance of about 55 percent is observed at the Curie temperature. 

Figure [7^ a) shows the magnetoresistance measured at 2 K up to an applied magnetic 
field of 3 kOe. It shows a maximum magnetoresistance of about 5.5 %. The change in 
magnetoresistance around 1 kOe is quite sharp and coincides with the SL phase indicated 
by magnetization measurements. At 2 K a small hysteresis is observed [upper inset in the 
Figure [7]^a)] and it appears in the same region where hysteresis is seen in the magnetization 
measurements [Fig. [2](a)]. The hysteresis persists up to 100 K (not shown). The lower 
inset shows the derivative of the resistivity with respect to the applied magnetic field which 
clearly shows the sharp change occurring at the applied field of 1 kOe. This behavior 
in the resistivity also seems consistent with the soliton model. As discussed above, as 
the applied magnetic field starts aligning the moments towards its direction, spin disorder 
scattering decreases thereby reducing the electrical resistivity. As the magnetic field causes 
the transition from the soliton phase to the ferromagnetic state, all the moments are aligned 
in the direction of the field and thus no further change in the resistivity at a particular 
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temperature. The low field magnetoresistance behavior is found to be qualitatively similar 
at all temperatures measured below 100 K. Figure [7](b) shows the normalized resistivity 
defined by Ph/ph=o as a function of magnetic field. Above 120 K the resistivity varies little 
with field, while below 120 K the resistivity behavior is similar to that observed at 2 K [cf. 
Fig. |7Fa)]. Interestingly, at Tq = 120 K p(H) has a strong field dependence but, unlike at 
lower temperatures, it does not saturate near 1 kOe. The origin of this behavior is unclear, 
and highlights the lack of detailed understanding of the complex behavior of helimagnetic 
materials near their critical temperatures. 
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FIG. 7: (Color online) Low field magnetoresistance of Cr 1 / 3 NbS2 measured in the aft-plane with 
magnetic field applied parallel to the plane, a) Magnetoresistance measured at T = 2 K. The 
upper inset shows the change in slope in the magnetoresistance in the vicinity of the field where 
the metamagnetic transition is observed in the magnetization measurements. The lower inset shows 
the derivative of the resistivity with respect of the field. It shows that the sharp change occurs at 
1 kOe. b) Normalized resistivity measured in the aft-plane as a function of magnetic field applied 
parallel to the plane at indicated temperatures. 
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The temperature dependent thermal conductivity k(T) of Cri/ 3 NbS2 is shown in figure 
(8^a). Upon heating, n increases up to about 65 K above which it decreases slightly up 
to 200 K. In general, the thermal conductivity of a metal is the sum of electronic and 
lattice terms. The electronic contribution to the thermal conductivity can be estimated 
by the Wiedemann- Franz Law^: K e p/T = L, where p is the electrical resistivity and L = 
2.45xl0~ 8 WfiK~ 2 is the Lorenz number. The inset in figure [8^a) shows the calculated 
electronic part of the thermal conductivity (n e ) for Cr^M^. This estimation suggests 
that the thermal conductivity is mainly due to phonons. 
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FIG. 8: (Color online) a) Temperature dependence of thermal conductivity measured in the ab- 
plane. Inset shows the electronic contribution to the thermal conductivity estimated from the 
Wiedemann-Franz law, and b) Temperature dependence of Seebeck coefficient measured in the 
a6-plane. Inset shows the magnetic field dependence of the Seebeck coefficient at T = 116 K with 
magnetic field applied parallel to the afr-plane 

Figure ^b) shows the temperature dependence of the Seebeck coefficient (S). At all 
temperatures measured, the Seebeck coefficient is positive. However, it shows a remarkable 
minimum slightly below the magnetic transition temperature. The Seebeck coefficient is very 
sensitive to the details of the band structure near the Fermi energy Ep. The anomalous 



13 



behavior observed near Tq may be related to the significant decrease in the density of states 
(DOS) around Ep in the magnetic state as compared to that in the non-magnetic state 
as presented below. Interestingly the Seebeck coefficient is found to have magnetic field 
dependence at 116 K. With the external magnetic field applied perpendicular to the c axis, 
S decreases up to about 20 kOe, above which it remains constant [Fig. fflb)] which suggests 
either a spin-entropy contribution to S or the effects of the magnetic field on the scattering 
of carriers. 



C. Heat capacity 

The temperature dependence of specific heat shows clear lambda anomaly in the vicinity 
of the magnetic ordering temperature at ambient field as shown in figure (9^a). The magnetic 
field dependence of the heat capacity anomaly is depicted in [Fig. [91(b)]. The lambda peak 
does not change significantly up to 800 Oe. However, it is found to decrease clearly at a 
field of 1500 Oe. The peak is suppressed with the further increase in the magnetic field and 
is not observed above 50 kOe. Application of a magnetic field increased the heat capacity 
above the magnetic transition temperature (120 K) and suppressed the heat capacity below 
it. A decrease in the lambda peak is found to occur in the ferromagnetic state when the 
moments are aligned in the direction of the magnetic field. This behavior is consistent with 
the results obtained in the magnetization and magnetoresistance measurements. 

The low temperature specific heat can be modeled well with electronic and phononic 
contributions. The inset of figure [9^a) shows a fit to C/T = 7 + /3T 2 , where, 7 is the 
Sommerfeld coefficient and /3 is the phononic heat capacity coefficient. The fit yielded 7 = 
0.40 ± 0.25 mJ mot 1 K~ 2 and (3 = 0.0026 ± 0.00047 mJ mot 1 K~ A . The electron density 
of states at the Fermi energy (Ep) and the Debye temperature (Tp>) can be estimated from 7 
and (3 using the relations N(Ep) = Sj/ir 2 )^ and Tp, = (5/3/12i?7r 4 ) _1//3 , respectively, where 
ke is the Boltzmann constant and R is the universal gas constant. The calculated values 
are N(Ep) = 3.40 states (eV-u.c.) -1 and Tp = 419 K. The experimentally determined DOS 
at the Fermi energy is comparable to the value [2.76 states (eV-u.c.) -1 ] obtained from the 



electronic structure calculations in the magnetic state presented below (Fig. 12). 
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FIG. 9: (Color online) a) Heat capacity of Cr^Nt^ as a function of temperature at zero applied 
magnetic field, b) Temperature dependence of the specific heat near the magnetic transition 
temperature measured at indicated magnetic fields applied perpendicular to the c axis. 



D. First principles calculations 

In order to understand theoretically the effect of the magnetic ordering we have per- 
formed first principles density functional theory calculations, using the generalized gradient 
approximation of Perdew, Burke and ErnzerhoP^ as implemented in the all electron code 
WIEN2k. 41 We have used the experimental hexagonal lattice constants and internal pa- 
rameters; no relaxation was performed. Calculations in both a non-magnetic state and a 
magnetic spin-polarized state were performed, using 120 /c-points in the irreducible Brillouin 
zone and LAP W sphere radii of 2.07 a , 2.33 ao and 2.34 ao for the S, Nb and Cr atoms 
respectively (here a is the Bohr radius = 0.529177 A). An RK max of 7, where R is the 
minimum atomic sphere radius and K the largest plane wave vector used in the expansion, 
was used. 

Regarding the spin-polarized calculations, a brief discussion is in order. As are previous 



15 



experimental evidence d 21 * 24 ^ that the ground state of Cr^Nh^ is a very long wavelength 
(~ 480 A) spiral with adjacent spins in the spiral very nearly parallel, such a magnetic 
structure, if one attempted to study it in a "brute-force" computational framework, would 
result in a calculation involving several hundred atoms and several thousand electrons, which 
would present a nearly intractable problem for standard first principles approaches. However, 
a simple observation allows one to perform a much simpler calculation which should be very 
close in energy and ground state properties to the actual ground state. The origin of the 
Dzyaloshinsky-Moriya interaction lies in the spin-orbit interaction] 42 * 43 ! The energy scale 
of the relativistic term believed to be responsible for the spiral state is much weaker (a 
quantitative estimate suggests between one and two orders of magnitude smaller) than the 
ordinary ferromagnetic Heisenberg nearest neighbor exchange energy. Thus it is highly likely 
that the actual energies and properties of the real ground state can be described as very 
nearly that of an ordinary spin-polarized ferromagnetic ground state, with the understanding 
that the "spin-up" and "spin-down" description of the properties presented below, such as 
band structure and density-of-states, are in essence only labels and that the actual physical 
state is a non-collinear spiral. It is then an excellent approximation to take the true ground 
state physical observable, such as the density-of-states, as the sum of the separate "spin- 
up" and "spin-down" contributions. Again, a quantitative estimate suggests that the effect 
of the DM interaction on the Kohn-Sham eigenvalues is less than 10 meV, which is much 
smaller than all other relevant energy scales and virtually invisible on the plots presented 
below. 

We begin with the calculated non-magnetic band structure and density-of-states, pre- 



sented in Figure 10 below. There are several bands crossing the Fermi level, and in addition 
flat bands virtually abutting the Fermi level from M to K and H to M. There is also a fairly 
complex structure around the T point. As might be expected this band structure leads to a 
high density-of-states at the Fermi level, presented in Figure [To|(b). The Fermi level sits near 
the middle of a region of fairly elevated DOS, with Fermi level DOS of approximately 15 
states/eV-unit cell and approximately 10 states/eV-unit cell Chromium. Given the two Cr 
in the unit cell and the exchange correlation value / of 0.38 eV for CrS2, one finds the Stoner 
criterion N I > 1 well satisfied, with N I on a per Cr basis having the value 1.9. The ma- 
jority, but not entirety, of the DOS character near Ep is Cr, suggesting some hybridization 
of the chemical bonds. 
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FIG. 10: (Color online) a) The calculated band structure b) and density of states of Cr 1 / 3 NbS2 
in the nonmagnetic state. Note the many Fermi level crossings and corresponding high density of 
states at Ep. 



Moving to the magnetic state calculations previously described, we find strong evidence 
for a magnetic ground state, with the magnetic state some 930 meV per formula unit in 



energy lower than the non-magnetic state. In Figure 11 we present the calculated band 
structure and density-of-states for this ground state; as outlined above we have simply 
summed the "spin-up" and "spin-down" contributions to the DOS and plotted both band 
structures on the same plot. The band structure demonstrates far fewer Fermi level crossings, 
and accordingly examining the density of states we find a substantial loss of spectral weight 
around Ep. It is instructive to directly compare the magnetic and non-magnetic DOS and 



this is presented in Figure 12 We note that, the significant changes to the DOS are confined 
to an energy range of ± 2 eV around Ep, despite the large energy gain of nearly an eV per 
formula unit. The DOS at Ep in the magnetic state is found to be 2.76 states (eV-u.c.) -1 
which is consistent with the experimental value extracted from the low temperature heat 
capacity presented above. 
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FIG. 11: (Color online) a) The calculated band structure and b) density of states of Cr 1 / 3 NbS2 in 
the assumed magnetic state. Most of the Fermi level crossings have disappeared and the density- 
of-states shows a substantial pseudogap at Ep. 
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FIG. 12: (Color online) The calculated density of states of Cr 1( / 3 NbS2 in the magnetic and non- 
magnetic states. 

V. SUMMARY AND CONCLUSIONS 

We have determined the temperature and magnetic field dependence of the magnetic, 
transport, and thermal properties of Cri/3NbS2 single crystals. Contours of constant mag- 
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netization are utilized to identify key features in the magnetic phase diagram near the 
transition region (Fig. [5]). The material is magnetically ordered below about 120 K. Below 
this temperature, increasing the magnetic field applied in the a6-plane results in a transition 
from helimagnetic to ferromagnetic order, with evidence from the previously reported soli- 
ton lattice phase at intermediate fields.^ Evidence for these field- and temperature- induced 
phase transitions is also seen in resistivity, Seebeck coefficient, and heat capacity results. 
The resistivity responds strongly to the magnetic ordering, decreasing upon cooling through 
the magnetic transitions by more than a factor of four. A sharp drop in the Seebeck coef- 
ficient is observed upon cooling near 120 K as well. Significant effects of applied magnetic 
field on these transport properties are also seen in the ordered state, especially near the 
ordering temperature. Observations include large magnetoresistiance (-55% at 140 kOe and 
120 K) and p(H) data that mimic the M(H) behavior below 120 K, but are anomalous at 120 
K. These findings suggest that the magnetic transitions are accompanied by large changes 
in the electronic structure, which is confirmed by DFT calculations. The calculated DOS in 
the magnetic state is strongly suppressed relative to the value from non-magnetic calcula- 
tions. Since the resistivity decreases upon entering the magnetically ordered state, changes 
in magnetic scattering of charge carriers are also likely to be important in determining the 
transport behavior. The large change in the resistivity suggests strong scattering by spin 
fluctuations above 120 K. Analysis of the heat capacity results indicate a small Sommerfeld 
coefficient of 0.4 mJmol^ 1 K^ 2 , in agreement with the magnetic DFT calculations, and a 
Debye temperature of 419 K. A lambda-like anomaly in the heat capacity near the magnetic 
transition is suppressed strongly at fields above about 1.5 kOe, when the material is in the 
ferromagnetic state. 
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